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Abstract 

The results of the first search for long-lived gluinos produced in 7 TeV pp collisions 
at the CERN Large Hadron Collider are presented. The search looks for evidence of 
long-lived particles that stop in the CMS detector and decay in the quiescent peri- 
ods between beam crossings. In a dataset with a peak instantaneous luminosity of 
1 x 10 32 cm _2 s _1 , an integrated luminosity of 10 pb -1 , and a search interval corre- 
sponding to 62 hours of LHC operation, no significant excess above background was 
observed. Limits at the 95% confidence level on gluino pair production over 13 orders 
of magnitude of gluino lifetime are set. For a mass difference nig — m^> > 100 GeV/ c 2 , 

and assuming BR(g — > gx\) = 100%, trig < 370 GeV/ c 2 are excluded for lifetimes from 
10 ^s to 1000 s. 
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Many extensions of the standard model predict the existence of new heavy quasi-stable par- 
ticles (J. Such particles are present in some supersymmetric models EHU, "hidden valley" 
scenarios 0, and grand-unified theories (GUTs), where the new particles decay through di- 
mension five or six operators suppressed by the GUT scale [6]. Long-lived particles are also 
a hallmark of split supersymmetry [7], where the gluino (g) decay is suppressed due to the 
large gluino-squark mass splitting, from which the theory gets its name. Of these possibilities, 
the Compact Muon Solenoid (CMS) experiment is most sensitive to models like split super- 
symmetry where production proceeds via the strong interaction resulting in relatively large 
cross sections at the Large Hadron Collider (LHC) |8l-[TT1|. For this reason, we have targeted the 
search described in this Letter at long-lived gluinos. Existing experimental constraints on the 
lifetime of such gluinos are weak fl2l H3| ; these gluinos may be stable on typical CMS exper- 
imental timescales. Lifetimes of 0(100-1000) seconds are especially interesting in cosmology 
since such decays would affect the primordial light element abundances, and could resolve 
the present discrepancy between the measured 6 Li and 7 Li abundances and those predicted by 
conventional big-bang nucleosynthesis H6llT4|. 

If long-lived gluinos were produced at the LHC, they would hadronize into gg, gqq, gqqq states, 
collectively known as "R-hadrons" some of which would be charged, while others would be 
neutral. Those that were charged would lose energy via ionization as they traverse the CMS 
detector. For slow R-hadrons, this energy loss would be sufficient to bring a significant frac- 
tion of the produced particles to rest inside the CMS detector volume ||T5). These "stopped" 
R-hadrons may decay seconds, days, or even weeks later, resulting in a jet-like energy deposit 
in the CMS calorimeter. These decays will be out-of-time with respect to LHC collisions and 
may well occur at times when there are no collisions in CMS. The observation of such decays, 
in what should be a "quiet" detector except for an occasional cosmic ray, would be an unam- 
biguous discovery of new physics. 

The CMS apparatus has an overall length of 22 m, a diameter of 15 m, and weighs 14000 
tons. The CMS coordinate system has the origin at the center of the detector. The z-axis points 
along the direction of the counterclockwise beam, with the transverse plane perpendicular to 
the beam; cp is the azimuthal angle in radians, 6 is the polar angle, and the pseudorapidity is 
rj = — In (tan [0/2]). The central feature of the CMS apparatus is a superconducting solenoid 
of 6 m internal diameter, providing a field of 3.8 T. Within the field volume are the silicon 
pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL) and the brass-scintillator 
hadronic calorimeter (HCAL). Muons are measured in gas-ionization detectors embedded in 
the steel return yoke. The HCAL, when combined with the ECAL, measures jets with a resolu- 
tion AE/ E 100 %y/ E [GeV] © 5 %. The calorimeter cells are grouped in projective towers, of 
granularity At] x Acp = 0.087 x 0.087 at central rapidities. In this analysis, jets are reconstructed 
using an iterative cone algorithm with R = \J Aip 1 + At] 2 = 0.5. The reconstructed jet energy E 
is defined as the scalar sum of the calorimeter tower energies inside the jet. The first level (LI) 
of the CMS trigger system, composed of custom hardware processors, uses information from 
the calorimeters and muon detectors to select (in less than 3 fis) the most interesting events. 
The High Level Trigger (HLT) processor farm further decreases the event rate to 300 Hz before 
data storage. A more detailed description of the CMS experiment can be found elsewhere |[T6]|. 

The 7 TeV center-of-mass pp collision data analyzed in this Letter were recorded by CMS be- 
tween April and October 2010. We divide these data into two samples: the first corresponds 
to 95 hours of trigger live-time during LHC "fills", in which the instantaneous luminosity was 
2 — 7 x 10 27 cm _2 s . We use this as a control sample to estimate the background rate. Because 
these data were recorded at relatively low instantaneous luminosity, there is negligible risk that 
a stopped-particle signal is present in this sample. The second sample, in which we search for 
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the presence of a stopped-particle signal, corresponds to 62 hours of trigger live-time during 
which data, corresponding to an integrated luminosity of 10 pb _1 , were recorded by CMS with 
a peak instantaneous luminosity of 1 x 10 32 cm _2 s _1 . In producing these data, the LHC was 
filled with up to 312 proton bunches per beam. 
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Figure 1: Probability for a produced R-hadron to stop anywhere inside the CMS detector for 
different gluino masses and -*/s = 7 TeV. The solid line shows the stopping probability for the 
cloud model with both electromagnetic (EM) and nuclear interactions (NI), the dashed line is 
for the EM only model, and the dotted line is for the "neutral R-baryon" model in which only 
R-mesons stop 07). Charge exchange reactions are considered NI. 



We employed a dedicated trigger to search for decays of particles at times when there are 
no collisions. Information from the beam position and timing (BPTX) monitors are used to 
identify gaps between the proton bunches that comprise the LHC beam. The BPTX monitors 
are positioned 175 m from the center of CMS on either side of the CMS interaction region and 
produce a signal when an LHC proton bunch passes the monitor. Even though the R-hadron 
decay does not produce a true jet, the resultant energy deposition is sufficiently jet-like that a 
jet trigger is reasonably efficient. We therefore require a jet trigger together with the condition 
that a coincidence of signals from both BPTX did not occur, ensuring that the trigger will not 
fire on jets produced from pp collisions. The LI trigger requires a jet with at least 10 GeV 
transverse energy. A 20 GeV threshold on jet energy is applied in the HLT At both LI and HLT, 
the pseudorapidity of the jet, |?/jet|/ is required to be less than 3.0. 

Additional selection criteria are applied during data analysis. Despite the BPTX veto in the 
trigger, several beam-related processes remain possible sources of background. In order to 
reject background events due to an unpaired proton bunch passing through CMS, events in 
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which either BPTX is over threshold are vetoed. Instrumental effects during trigger generation, 
and features of the LHC beam such as lower intensity "satellite" bunches that accompany the 
colliding protons, can cause triggers in some of the 25 ns intervals (BX's) which precede or 
follow the one in which the intended proton collisions occur. We therefore reject any event oc- 
curring up to two BX's before, or one BX after, the BX in which collisions are expected. To reject 
beam-halo muon events, which may not be synchronous with proton collisions, we employ a 
loose beam-halo veto using the cathode strip chambers (CSC) in the endcap muon system. The 
algorithm rejects events in which a beam-halo trigger was recorded, or a track segment was 
reconstructed in the CSC system with timing consistent with beam-halo, or a muon track was 
reconstructed with beam-halo-like kinematics. Finally, to ensure that no out-of-time pp colli- 
sion events due to satellite bunches contaminate the search sample, events with one or more 
reconstructed primary vertices are rejected. 

A small fraction of cosmic rays traversing the CMS detector deposit significant energy in the 
calorimeters. To reduce this background, events that contain reconstructed muons are vetoed. 
Once beam-related backgrounds and cosmic rays are removed, the remaining source of back- 
ground is instrumental noise. Standard calorimeter cleaning and noise rejection criteria [T181 - I20I 
are applied. We restrict our search to jets in the central HCAL; we require that the most ener- 
getic jet in the event has \r/^ et \ < 1.3. To suppress noise fluctuations and energy deposits from 
cosmic rays, a jet with reconstructed energy above 50 GeV is required. To remove events where 
a single HCAL channel has misfired, events with more than 90% of the energy deposited in 
three or fewer calorimeter towers are vetoed. We also require that the leading jet has at least 
60% of its energy contained in fewer than 6 towers. To suppress noise from hybrid photodiode 
discharge fl8 | , events with 5 or more of the leading towers at the same azimuthal angle, or 
where more than 95% of the jet energy is contained within towers at the same azimuthal angle, 
are rejected. 

The HCAL electronics have a well-defined time response to charge deposits generated by 
showering particles. Analog signal pulses produced by these electronics are sampled at 40 MHz, 
in synchrony with the LHC clock. These pulses are readout over ten BX samples centered 
around the pulse maximum. A physical pulse has some notable properties which we use to dis- 
tinguish it from noise pulses. There is a clear peak in the signal pulse shape (BX pea ]J, significant 
energy in one bunch crossing before the peak (BX pea k_i), and an exponential decay for several 
BX's following the peak. We use the ratios R\ = BX pea ^ + i / BX pea ^ and R2 = £>X pea k + 2/ ^^peak+i 
to characterize the exponential decay, requiring Ri > 0.15 and 0.10 < R2 < 0.50. Since a physi- 
cal pulse spans only four time samples we are able to reject noise events based on the presence 
of energy in previous or successive BX's. We remove events with more than 10% of the energy 
of the pulse outside of the central four BX's. Energy deposits from physical particles tend to 
have a large fraction of the pulse energy in the peak BX. Noise does not produce this pulse 
shape; noise pulses tend to be spread across many BX's or localized in one BX. We require 
the ratio of the peak energy to the total energy to be between 0.4 and 0.7. The requirement 
that the pulse shape satisfy all the preceding criteria rejects 50% of the remaining events in the 
background sample while preserving 93% of simulated signal events. 

We have developed a custom, factorized simulation of gluino production, stopping, and decay 
to investigate the experimental signature of this atypical signal. First, we generate qq — > gg 
and gg -4 gg events at y/s = 7 TeV using Pythia |2l). The lifetime of the gluino is set such 
that it is stable. Gluino masses ma = 150 to 500 GeV/c 2 are studied. Pythia hadronizes the 
produced gluino into R-hadrons. A modified Geant4 [22J that implements a "cloud model" of 
heavy stable colored interactions with matter [23] is used to simulate the interaction of these R- 
hadrons with the CMS detector and to record the location at which those R-hadrons that do not 
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exit the detector come to rest. Figure [T] presents the stopping probability as a function of gluino 
mass obtained from our simulation for this model and two alternative models of R-hadronic 
interactions with matter. Next, we again use Pythia to produce an R-hadron at rest which 
we translate from the nominal vertex position to the recorded stopping location and decay the 
constituent gluino instantaneously via g — > g X.\- Finally we use a specialized Monte Carlo 
simulation to determine how often the stopped gluino decay would occur during a triggerable 
beam gap. Further details of this simulation are described elsewhere [24 j. 

The efficiency with which triggered events pass all selection criteria is estimated from the sim- 
ulation to be 54% for a representative gluino decay signal (nig = 300 GeV/c 2 and m^o = 

200 GeV/c 2 ). This point was chosen to be above existing limits and within the reach of CMS. 
The equivalent efficiency with respect to all stopped particles is 17% since a significant num- 
ber of R-hadrons stop in uninstrumented regions of the CMS detector where their subsequent 
decay would not be observable. For any new physics model that predicts events with suf- 
ficient visible energy trig — m^o > 100 GeV/c 2 , this efficiency does not change significantly. 
We measure the background rate in the control sample after all but one of the selection cri- 
teria are applied, R£cmtror ^ e a ^ so measure the background rate in the control sample af- 
ter all selection criteria are applied, R^ ntrol - To obtain an estimate of the background rate in 
the search sample after all selection criteria, we again measure the rate after omitting one se- 
lection criterion and multiply it by the ratio of the rates obtained from the control sample, 
R s^arch = ^^(^ntrol/^nD- This procedure is performed twice, each time omitting one 
of the most powerful background rejection criteria such that R N_1 3> R N ; we take the mean of 
both determinations as the final background rate estimate. 

We estimate the systematic uncertainty on the background to be 23% from the observed vari- 
ation of the control sample rate R^troi during the time period in which the data were taken. 
There is also a potential systematic uncertainty due to the accuracy with which the energy de- 
position of our jet-like signal is simulated. From proton/ pion test-beam data and studies of the 
energy deposited in HCAL by incident cosmic rays, we estimate this introduces a 7% uncer- 
tainty on the acceptance. The systematic uncertainty due to trigger efficiency is negligible since 
the data analyzed are well above the turn-on region. Similarly, the systematic uncertainty due 
to reconstruction efficiency is negligible since we restrict our search to nig — m^o > 100 GeV/ c 2 
wherein we are fully efficient. Finally, there is an 11% uncertainty on the luminosity measure- 
ment [25]. Limits on a particular model (e.g., gluinos in split supersymmetry) introduce more 
substantial systematic uncertainties, since the signal yield is sensitive to the stopping probabil- 
ity. The stopping probability varies greatly depending on the model of R-hadronic interactions 
used in the simulation. 

After the selection criteria described in the preceding paragraphs are applied, we perform a 
counting experiment and a time-profile analysis on the remaining data. For the counting ex- 
periment, we consider gluino lifetime hypotheses from 75 ns to 10 6 seconds, where we have 
chosen the upper limit of the search to be the longest lifetime for which we can still expect 
to observe at least one event. For lifetime hypotheses shorter than one LHC orbit (89 jia), we 
search within a time window following each filled bunch crossing. This time window is equal 
to 1.256 x Tg for optimal sensitivity to each hypothesized gluino lifetime Xs. In addition to the 
lifetimes required to map the general features of the exclusion limit, we include two lifetimes 
for each observed event: the largest lifetime hypothesis for which the event lies outside the time 
window, and the smallest lifetime hypothesis for which the event is contained within the time 
window. For lifetime hypotheses longer than one LHC fill, we do not consider the possibility 
that any observed events may have come from gluinos produced in a previous fill. 
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In the search sample, we do not observe a significant excess above expected background for any 
lifetime hypothesis. The results of this counting experiment for selected lifetime hypotheses are 
presented in Table [l] In the absence of any discernible signal, we proceed to set 95% confidence 
level (C.L.) limits over 13 orders of magnitude in gluino lifetime using a hybrid CLg method [26 1 
inspired by Ref. |I27| . In Fig.|2]we show the 95% C.L. limit on o~(pp — > gg) x BR(g — > gx®) for 



a mass difference trig 



> 100 GeV/c 2 . The error bands include statistical and systematic 
uncertainties. With the horizontal line in Fig. [2] we show a recent NLO+NLL calculation of the 
cross section at \/s = 7 TeV for nis = 300 GeV/c 2 from the authors of Ref. [11]. To illustrate 
the effect of the stopping probability uncertainty, we present three different 95% C.L. limits on 
o~(pp — > gg) x BR(g — > gXi) i n which the three different R-hadron models are used. Assuming 
the cloud model for the interaction of R-hadrons with matter, and assuming BR(g — > gx®) = 
100%, we are able to exclude lifetimes from 75 ns to 3 x 10 5 s for nig = 300 GeV/c 2 with the 
counting experiment. Finally, we present the result as a function of the gluino mass in Fig. [3] 
Under the same assumptions as for the cross section limit, we exclude nig < 370 GeV/c 2 for 
lifetimes between 10 f/s and 1000 s. If we assume the EM only model for R-hadronic interactions 
with matter in order to compare with what was done in Ref. [12], this exclusion becomes nig < 
302 GeV/c 2 . 

Table 1: Results of counting experiments for selected values of Tg. Entries between 1 x 10~ 5 
and 1 x 10 6 s are identical and are suppressed from the table. 



Lifetime [s] 


Expected Background (± stat. ± syst.) 


Observed 


1 x 10" y 


0.8 ±0.2 ±0.2 


2 


1 x 10~ 6 


1.9 ±0.4 ±0.5 


3 


1 x 10~ 5 


4.9 ±1.0 ±1.3 


5 


1 x 10 6 


4.9 ±1.0 ±1.3 
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We also perform a time-profile analysis. Whereas, for short lifetimes, a signal from a stopped 
gluino decay is correlated in time with the collisions, backgrounds are flat in time. Since the 
signal and background have very different time profiles, it is possible to extract both their con- 
tributions by analyzing the distribution of the observed events in time. We assume all colliding 
bunches in an orbit have equal individual instantaneous luminosity. We build a probability 
density function (PDF) for the gluino decay signal as a function of time for a given gluino life- 
time hypothesis and the actual times of LHC beam crossings as recorded in our data. Figure [4] 
shows an example of such a PDF for a gluino lifetime of 1 }is; the in-orbit positions of 2 ob- 
served events in the subset of our data that were recorded during an LHC fill with 140 colliding 
bunches are overlaid. We limit the range of lifetime hypotheses considered for this time-profile 
analysis to 75 ns to 100 jis such that the gluino lifetime is not much longer than the orbit period. 
For each lifetime hypothesis we build a corresponding signal time profile, fit the signal plus 
background contribution to the data, and extract a 95% C.L. upper limit on the possible signal 
contribution. The obtained results are plotted as a dotted line in Fig. [2] This temporal analysis 
relies only on the flatness of the background shape; it does not have the counting experiment's 
systematic uncertainty on the background normalization. Consequently, its dominant system- 
atic uncertainty is the 11% uncertainty on the luminosity measurement. For a mass difference 
nig — m^o > 100 GeV/c 2 , assuming BR(g — > gx®) = 100%, we are able to exclude nig < 382 

GeV/ c 2 at the 95% C.L. for a lifetime of 10 ]is with the time-profile analysis. 

We have presented the results of the first search for long-lived gluinos produced in 7 TeV pp 
collisions at the LHC. We looked for the subsequent decay of those gluinos that would have 



6 



_Q 
Q. 



T 

DC 
CD 

X 

lS 

T 

Q. 





1 — 1 — 1 — 1 — 1 

■ J L dt = 1 pb" 1 


i i — i — i — i — i — i — i 

95% C.L. Limits: 

Expected: Counting Exp. 


10 4 


r C x = 1 x 10 32 cm- 2 s- 1 


Expected ±1a: Counting Exp. 
txpectea ±^o. L-ounting txp. 




; Vs = 7 TeV 


Obs.: Counting Exp. 

Obs.: Counting Exp. (Neutral R-Baryon) ■ 


10 3 


m~ - nrLo = 1 00 GeV/c 2 
_ g z 

"*>. s 


Obs.: Counting Exp. (EM only) 

Obs.: Timing Profile /- 

/ 

/ 

/ / 


10 2 




/ . 
/ i 

NLO+NLL (m = 300 GeV/c 2 ) / / 
9 /: 








w J : 






10 


■i i i i i 


i i i i i i i i " 



10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 1 10 10 2 10 3 10 4 10 5 10 6 

Figure 2: Expected and observed 95% C.L. limits on gluino pair production cross section times 
branching fraction using the "cloud model" of R-hadron interactions as a function of gluino 
lifetime from both the counting experiment and the time-profile analysis. Observed 95% C.L. 
limits on the gluino cross section for alternative R-hadron interaction models are also pre- 
sented. The NLO+NLL calculation is for nig = 300 GeV/ c 2 from a private communication 
with the authors of Ref. IfTTl. 



stopped in the CMS detector during time intervals where there were no pp collisions. In par- 
ticular, we searched for decays during gaps in the LHC beam structure. We recorded such 
decays with dedicated calorimeter triggers. In a dataset with a peak instantaneous luminosity 
of 1 x 10 32 cm _2 s _1 / an integrated luminosity of 10 pb -1 , and a search interval corresponding 
to 62 hours of LHC operation, no significant excess above background was observed. Limits 
at the 95% C.L. on gluino pair production over 13 orders of magnitude of gluino lifetime are 
set. For a mass difference nig — m^o > 100 GeV/c 2 , assuming BR(g — > gx\) = 100%, we ex- 
clude nig < 370 GeV/c 2 for lifetimes from 10 to 1000 s with a counting experiment. Under 
the same assumptions, we are able to further exclude nig < 382 GeV/c 2 at the 95% C.L. for a 
lifetime of 10 jis with a time-profile analysis. These results extend existing limits from the DO 
Collaboration |Tl2| on both gluino lifetime and gluino mass. These limits are the most restrictive 
to date. 

We wish to congratulate our colleagues in the CERN accelerator departments for the excellent 
performance of the LHC machine. We thank the technical and administrative staff at CERN and 
other CMS institutes, and acknowledge support from: FMSR (Austria); FNRS and FWO (Bel- 
gium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and 
NSFC (China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); Academy of Sci- 
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Figure 3: 95% C.L. limits on gluino pair production cross section times branching fraction as 
a function of gluino mass assuming the "cloud model" of R-hadron interactions (solid line) 
and EM interactions only (dot-dashed line). The nig — m^o mass difference is maintained at 100 

GeV/ c 2 ; results are only presented for m^o > 50 GeV/ c 2 . The NLO+NLL calculation is from a 



private communication with the authors of Ref. ffTT|. The lifetimes chosen are those for which 
the counting experiment and time-profile analysis are most sensitive 
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